Possible Stratification Mechanism in Granular Mixtures 
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We propose a mechanism to explain what occurs when a mixture of grains of different sizes 
and different shapes (i.e. different repose angles ) is poured into a quasi-two-dimensional cell. 
Specifically, we develop a model that displays spontaneous stratification of the large and small 
grains in alternating layers. We find that the key requirement for stratification is a difference in 
the repose angles of the two pure species, a prediction confirmed by experimental findings. We also 
identify a kink mechanism that appears to describe essential aspects of the dynamics of stratification. 

PACS Numbers: 05.40+j, 46.10+z, 64.75+g 
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Granular materials exhibit many unusual properties 
, such as size segregation, when exposed to external vi- 
brations or rotations |2) . Size segregation is also observed 
when a mixture of grains of different size is poured onto a 
pile || , the large grains are found preferentially near the 
bottom of the pile, while the small grains are found near 
the top. Recently, it was found B that when a mixture 
of grains of different sizes and shapes is poured between 
two vertical slabs separated by « 5 mm, there appears 
a spontaneous stratification, with alternating layers of 
small and large grains parallel to the surface of the sand- 
pile. Additionally, there is an overall tendency for the 
large and small grains to spontaneous segregate in differ- 
ent regions of the cell j|[IJ. 

Very recently, Boutreux and de Gennes (BdG) || 
treated theoretically the case of granular flow made of 
two species. They based their work on a set of coupled 
convection equations to govern the flow of rolling grains 
and their interaction with the sandpile, introduced earlier 
by Bouchaud et al. in the case of a single-species sandpile 
H . BdG reproduced the phenomenon of segregation, but 
an understanding of stratification is lacking. 

Here we seek to understand segregation and stratifica- 
tion in the conditions of Q, where the two species have 
different size and different shape. We first introduce a 
discrete model to give a clear picture of the phenomenol- 
ogy, and then develop a continuum approach. In agree- 
ment with the experimental findings pj , we find that seg- 
regation is related to the difference of size of the grains, 
and stratification, to the difference in repose angles of 
the two pure species. 

In the discrete model, the sandpile is built on a lattice, 
where the grains have the same horizontal size as the lat- 
tice spacing and two different heights, Hi and H2 > H\. 
Each grain belongs to one of two phases: a static phase 
(if the grain is part of the solid sandpile) and a rolling 
phase (if the grain is not part of the sandpile but rolls 
downward with a constant drift velocity) (| . The local 
slope Si = hi — hi + i of the static grains is the variable 
controlling the dynamics of the rolling grains, where hi 
is the height of the sandpile at column i. 

At each time step, we deposit at the top of column 1 
of the pile Ni small grains plus N2 large grains; these 



grains belong to the rolling phase. One rolling grain per 
column of each species interacts with the surface at each 
time step, and can be converted from the rolling phase 
to the static phase. The remaining rolling grains are con- 
verted downward with unit drift velocity — i.e., they all 
move to the next column at each time step. 

The dynamics of each rolling grain interacting with 
the sandpile surface is governed by its repose angle (the 
maximum angle below which a rolling grain is converted 
into a static grain) H|). We note that the repose angle 
depends on the local composition on the surface, so we 
define 9 a /3 as the repose angle of a rolling grain of type 
a on a surface with local grains of type f3. We choose 
$2i < On to take into account that large grains roll more 
easily on top of small grains than small grains roll on top 
of large grains (since the surface "looks" smoother for 
large grains rolling on top of small grains, see Fig. |l]a). 
The repose angles of pure species 9 aa lie between #21 and 

#12- 

The stratification experiments use a mixture of 
grains of different shapes (smaller "less faceted" grains 
and larger "more faceted" grains). The repose angle of 
the smaller pure species is then smaller than the repose 
angle of the large pure species — i.e., 9\\ < 622- To mimic 
the experimental conditions for stratification [0], we set 
#21 < &ii < #22 < #12- We notice that the condition 
#21 < 612 is a consequence of the different size of the 
species, while the condition 6n ^ 6*22 is achieved, e.g., 
for mixtures of grains with different shapes. 

At each time step, the rolling grain interacting with 
the sandpile surface at coordinate i will either stop — by 
being converted into a static grain — if the local slope of 
the surface hi — /li+i < s a p = tan 9 a p, or will continue to 
roll (together with the remaining rolling grains) to col- 
umn i + 1 if hi — hi + \ > s a p. We iterate this algorithm 
to form a large sandpile of typically 10 5 grains. 

Figure [IJd shows the resulting morphology. The strat- 
ification is qualitatively the same as that found experi- 
mentally Q , not only in regard to the statics of the sand- 
pile (seen in Fig. [l]b), but also in regard to the dynamics. 
After a pair of static layers is formed with a layer of large 
grains on top of a layer of small grains, the angle of the 
sandpile is close to #22- Since the surface of the sand- 
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pile is made of large grains and #22 < #12, a thin layer of 
small grains is trapped on top of the layer of large grains. 
These small grains smooth the surface without changing 
significantly the sandpile angle, and allow rolling small 
grains to go further down (since On < 622)- The large 
grains are rolling down on this thin layer of small grains 
without being captured (821 < #22), and are the first to 
reach the bottom of the sandpile, giving rise to segre- 
gation. When the flow reaches the base of the pile, the 
grains develop a profile which displays a "kink" where 
the grains are stopped: the small grains are stopped first 
since 6*21 < 6*12, and a pair of layers begins to form, with 
the small grains underneath the large grains. The kink 
moves upward (in the opposite direction to the flow of 
grains) until it reaches the top of the pile and a complete 
pair of layers has been formed. 



J 12 > J 21 
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FIG. 1. (a) The slopes of repose s a /3 depend on the 
composition of grains at the surface of the pile, and are 
chosen according of the four possible interactions between 
small and large grains. The slopes of repose satisfy 
S21 < sii < S22 < S12 (see discussion in text), (b) Result ob- 
tained with the discrete model (for Hi = 1, H2 = 2, S21 = 2, 
sia = 6, S22 = 7, S12 = 10, Ni = 20, and N 2 = 10). 



If, on the other hand, we consider #22 < #11 in the 
model (corresponding to a mixture of smaller "more 
faceted" grains, and larger "less faceted" grains) we find 
segregation but no stratification. Thus, the control pa- 
rameter for stratification appears to be the difference in 
the repose angle of the pure species. 



As seen in Fig. |]b, before the layers appear there is 
an initial regime where only segregation is found. At 
the onset of the instability leading to stratification, a 
few large grains are captured on top of the region of 
small grains near the center of the pile where the an- 
gle of the pile is 9 ~ On- The repose angle for large 
grains is now #22- Thus, if 9 ~ 9n < 6*22, more large 
grains can be trapped (since the angle of the surface is 
smaller than the repose angle), leading to the first sub- 
layer of large grains and then to stratification. On the 
other hand, if ~ 9n > 6> 2 2, no more large grains are 
captured, the fluctuation disappears, and the segregation 
profile remains stable. Thus this picture suggests that, 
in agreement with Q , the segregation profile observed in 
the initial regime is "stable" so long as 6*22 < #11, and 
unstable (evolving to stratification for large enough sys- 
tems) when On < 6^2 ■ 

To offer further insight into the physical mechanism 
of stratification, we now develop a continuum approach 
in the spirit of Refs. The variables are the two 

thicknesses of rolling grains R a (x,t), with a = 1,2 re- 
spectively for small and large grains, the height of the 
sandpile h(x,t), and the volume fraction of grains of /? 
type in the static phase 4>p(x, t). The equations of motion 
are g§ 

dR a dR a 

1 a (Taj 



dt 
dh 

~dt 



dx 



(ib) 



Here v a is the downhill convection velocity of species 
a, and r a describes the interaction of the rolling grains 
with the surface — i.e., how rolling grains are stopped and 
become part of the sandpile (capture), and how grains 
of the sandpile can enter the rolling phase (amplifica- 
tion) . The concentrations are given by <f>i + <f>2 = 1 5 and 
cj> a = -T a /(dh/dt). 

As in the discrete model, we focus on the dependence 
of the repose angle on the composition of the surface 
<j)p{x,t). The repose angle 9 a of each type of rolling 
grain is now a continuous function of the composition 
of the surface 9 a = 9 a (4>p) (see Fig. ||a). The repose an- 
gle 9 a f3 defined for the discrete model is now 6 a ((pp) with 
<f)p = 1. We propose that T a = T a (R a , <pp, 9\ oc ) obeys 
the relation 



loc{9\oc 9 a 

la <t>a{9\oc - 9a(<fip)) Ra 



if 0i oc < 6 a (<f>p) 
if 9 loc > 9 a ((j)p) 



(2) 



where 9\ oc (x, t) = —dh/dx is the local surface angle. The 
parameter "f a represents the effectiveness of the interac- 
tion: v a jia ~ d a (where d a is the linear size of the grain) 
is the length scale on which a rolling grain will interact 
significantly with the surface when 9\ oc is slightly differ- 
ent from the repose angle 0. 

The interaction term T a includes two types of process 
||. (a) Capture: as in the discrete model, rolling grains 
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are captured if 8\ oc < a (<Pp)- The capture is propor- 
tional to the number of grains interacting with the sur- 
face, and so is proportional to R a . (b) Amplification: if 
#ioc > 0a(<t>/3), then some static grains of the sandpile are 
converted into rolling grains. This conversion of a type 
grains is proportional to their surface concentration </>„, 
and to the number R a of rolling grains acting in ampli- 
fication. 

We next solve Eqs. (0)-§) numerically fie] ]. The re- 
sults, shown in Fig. |^c, are qualitatively similar to the 
discrete model [jUI. We find stratification whenever the 
repose angle has the qualitative behavior shown in Fig. 
||a, the key requirement being 622 > #11 ■ We also find a 
"kink," corresponding to the growth of the new pair of 
static layers, with a well-defined steady-state profile and 
upward velocity. 



%jj = 6i((f>2) — ^2(^2) is independent of the concentration 
$2, and we set v\ — V2 = v, and 71 = 72 = 7 (Fig. ||b). 
We seek a steady-state solution, where the profiles of the 
sandpile and of the rolling grains are conserved in time. 
Thus, stratification cannot be observed for this solution. 
The conservation of the grains gives dh/dt = vR°/L, 
and we impose dR a /dt = 0, with boundary conditions 
R a (0) =R% = R°/2 and R a (L) = 0. 

The steady-state solution of (|l|)-(§) shows almost to- 
tal segregation. At the upper part of the pile, for < 
x < x mi with x rn = L/2 — v/l^ip), only small grains are 
present (tp^x) = 1, and 2 ( a; ) = 0), and the profiles are 
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FIG. 2. (a) Dependence of the repose angle for the two 
types of rolling grains on the concentration of the surface of 
large grains (j>2- An essential ingredient to obtain stratifica- 
tion is that 822 > On- For the numerical integration, we use 
the linear interpolation between <j)2 ~ and 02 = 1 as plotted 
here, (b) Picture with the different quantities appearing in the 
text. The dash-circled zone is the "kink." (c) Resulting mor- 
phology of the numerical integration of the continuum equa- 
tions. The parameters used are tan(9n = 1, tan O22 = 1-1, 
tan 6<i2 = 1.4, tan 6*21 = 0.7 , 71 = 72 = 0.8, and Vi = V2 = 1. 



To find the conditions under which stratification oc- 
curs, we first calculate the steady-state solution of the 
full set of Eqs. (|l|)-(^), and then study its stability 
under perturbations. To describe the experimental sit- 
uation, we consider a 2-D cell with vertical walls at 
x = and x = L B- We assume that the difference 



hoc(x) - On 



-v/l 



L/2 



(3b) 



At the lower part of the pile (x m < x < L), we find 
that, after a small region of size of the order of v/^ip), 
mainly large grains are present, and the profiles are 



4>i(x) — exp 
Ri{x) 



2v Mx)R(x) 



hoc{x) - O22 = -m<l>y{x) 



l(L-x)' 



(4a) 
(4b) 
(4c) 



Here R = i?i + R 2 = R°(l - x/L), and m = 6» 22 - O21 = 
612 — On- 

To analyze the stability of the steady-state solution (||) 
and (^) for the different phenomenologic parameters, we 
impose the steady-state solution as the initial condition, 
and then we look numerically for the stability of the pro- 
file under perturbations. For 6>n > 822, the steady state 
solution is stable: in this case, only segregation is ob- 
served, and the sandpile conserves in time the profiles 
(||) and (||). For On < 622, the steady-state solution is 
unstable (evolving to stratification), just as in p(. 

To gain insight about the "kink" mechanism, we look 
for a possible steady-state solution for the shape of the 
kink assuming that (i) far below and above the kink, the 
sandpile has a constant angle 6q\ (ii) the lowest part of 
the kink is made only of small grains, so that large rolling 
grains are not captured, and the top part is made only 
of large grains. 

To suppose the existence of a stationary solution for 
the kink implies that R\(x,t) and f(x,t) = h(x,t) + 8qx 
are functions only of u = x + Vjt, where is the uphill 
speed of the kink. For the lowest part of the kink, as only 
small grains are captured (</>i(u) = 1, i?2(u) = R°/2), 
Eqs. (|l|) reduce to equations for Ri(u) and f(u). We 
obtain the shape of the low part of the kink: for u < 0, 
f(u) = and for u > 0, f(u) obeys 
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•log 



(5) 



where <5i = 8q — 9\ i , and u> = u-j- / (u + u-f ) . Then the lower 
layer of the kink is characterized by a linear dependence 
f(u) oc u§i for u <C R /(2w5i), plus logarithmic correc- 
tions near the boundary with the upper layer of large 
grains. This solution is no longer valid when the angle 
of the surface reaches 6*21, and the large grains start to 
be captured. We note that this stationary solution exists 
only when <5i > 0. 

The solution of the equations for the highest part of the 
kink where only large grains are present can be obtained 
in the same way and is 



/(«) 



w 



1 — exp 



(6) 



where 82 = #0 — O22 ■ We then find that the shape of the 
upper part of the kink is exponential, and exists only for 
82 < 0. 

Thus we see that the existence of the stationary solu- 
tion for the kink implies that 9n < 9q < 622- the sandpile 
is built on an angle intermediate between the two repose 
angles of the pure species, and the repose angle of the 
small grains must be smaller than the repose angle of the 
large grains — in agreement with experiments and 
the stability analysis performed above. 

The layer thickness A is R°/w (see Eq. ^J), which is a 
consequence of the conservation law stating that all the 
rolling grains are stopped at the kink [Q. Furthermore, 
Eq. (ra) implies jivf/v^ ~ 1. For / ~ R°/w, this gives 
Uf oc 7-R , so that, for R° 0, we obtain 



A cx c u/7 + R°. 



(7) 



where c is a numerical constant that does not depend on 
v, 7, or R° . This relation, which we verify numerically, 
is relevant since v/7 and R° are both of the order of the 
diameter of the grains. 

The typical size of the initial regime of segregation, L x , 
observed prior to stratification when 8 = O22 — On > 
(Fig. ||c), can be calculated as follows. The condition for 
the appearance of a first layer of large grains on top of 
the region of small grains near the center of the pile is 
that capture of large grains must be larger than capture 
of small grains, i.e. |T2| > |Ti | , where (Til = jmRi, and 
T2 1 = 7<5i?2- Assuming that the solution (3|) is valid for 
the initial segregation regime, we can evaluate R%, and 
i?2 at x = x m . We obtain 



L r 



v m R 



(8) 



and verify (ph numerically. 

In conclusion, we develop a mechanism to explain the 
observed stratification H. This mechanism is related to 



the dependence of the local repose angle on the local sur- 
face composition. We find that stratification occurs only 
when the repose angle of the large grains is larger than 
the repose angle of the small grains (6*22 > #11, corre- 
sponding to large grains rougher than small grains). The 
model describes the static picture of the sandpile of Q| 
with alternating layers made of small and large grains, 
and also reproduces the dynamics, where the layers are 
built through a "kink" mechanism. When 622 < &11, 
the model predicts almost complete segregation, but not 
stratification. These results are in agreement with exper- 
iment [Q. 
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